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1. Introduction 



3. Event selection 



Positron-electron annihilation produces 
hadronic final states with an amplitude pro- 
portional to 

A ~ -j»J», (1) 
s 

where e is the charge of the electron, s is the 
square of the center-of-mass energy, is the 
e + e~ current, and is the hadronic current 
for the final state. The object of many experi- 
ments is to measure the matrix elements of J M . 
In e + e _ — > pp, a pair of spin-1/2 baryons with 
internal structure are produced, and the current 
contains two independent form factors, electric 
and magnetic, G^g 2 ) and Gm(? 2 ) 00- 

Understanding nucluon structure is one of the 
central problems of hadronic physics. In the 
time-like region, two processes, e + e~ — > pp and 
pp — ► e + e~, are used to measure the proton 
form factors Gs(g 2 ) and Gj./fe 2 ) as functions of 
the four- momentum transfer q 2 . Data samples 
in previous experiments are limited 0-0 • For 
6 < s < 8 GeV 2 , there is no data up till now. 

In this letter, we use the data from the up- 
graded Beijing Spectrometer (BESII) at the Bei- 
jing Electron-Positron Collider (BEPC) cover- 
ing the center-of-mass energy of 2.0-3.0 GeV in 
1999 8 and the data at 2.2, 2.6, and 3.07 GeV in 
2004 to measure the cross section of e+e~ — ► pp, 
and also determine the proton form factor in this 
energy range. 

2. BES detector 

BES is a conventional solenoidal magnet detec- 
tor that is described in detail in Ref. ;9^ ; BESII is 
the upgraded version of the detector [H]. A 12- 
layer vertex chamber (VC) surrounding the beam 
pipe provides trigger and coordinate information. 
A forty-layer main drift chamber (MDC), lo- 
cated radially outside the VC, provides trajectory 
and energy loss (dE/dx) information for charged 
tracks over 85% of the total solid angle. The mo- 
mentum resolution is a p jp = 0.017^/1 + p 2 (p 
in GeV/c), and the dE/dx resolution for hadron 
tracks is ~ 8%. An array of 48 scintillation coun- 
ters surrounding the MDC measures the time-of- 
fllght (TOF) of charged tracks with a resolution 
of ~ 200 ps for hadrons. Radially outside the 
TOF system is a 12 radiation length, lead-gas 

Kftrrnl QVinwpr r>mmtr>r fRSIf^ TlniQ mnflQnrr>Q tlnr> 



To select e + e — * pp, the following criteria are 
used: 

1. There must be two oppositely charged 
tracks in the MDC. Each track should have 
a good helix fit in the polar angle range 
| cos0| < 0.8 in the MDC, and the point of 
closest approach of the tracks to the beam 
axis should be within 2 cm in the radial di- 
rection and within 15 cm of the interaction 
point longitudinally. 

2. Tracks should be back-to-back. Fig. ^ 
shows the acollinearity (Acol) distributions 
for Monte Carlo (MC) simulated events at 
<Js = 2.0, 2.2, 2.6, and 3.0 GeV. Because 
of energy loss, the Acol distribution for low 
energy charged particles is somewhat differ- 
ent from higher energy ones. For 2.0 GeV 
data, we require that Acol is less than 10°, 
and for other energy points, less than 3°. 

3. To remove Bhabha events, E p /p p < 0.6 is 
required, where E p is the deposited energy 
in the BSC and p p is the momentum of the 
candidate proton. Fig. shows E p /p p dis- 
tributions for protons and positrons for MC 
e + e~ — * pp events and Bhabha events at y/s 
— 2.2, 2.6, and 3.0 GeV. This requirement 
removes most Bhabha background. 
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Figure 1. Acol distributions for MC e + e — > pp 
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Figure 2. Distributions of the ratio of the de- 
posited energy in the BSC and the momentum 
for final state particles of MC e + e~ — > pp events 
and Bhabha events at y/s = 2.2, 2.6, and 3.0 GeV. 
Left plots are protons, right are positrons. 



5. Lastly, dE/dx information is used to iden- 
tify pp pairs at y/s = 2.0 and 2.2 GeV, and 
TOF information is used at all other energy 
points JI] . Each charged track must satisfy 
Probp > 0.01, where Prob p is the particle 
identification confidence level for the proton 
or antiproton hypothesis determined cither 
from dE/dx or TOF information. 

The numbers of events passing the selection cri- 
teria are listed in Table 2] 

4. Luminosity 

The integrated luminosity is determined from 
large-angle Bhabha events using Q3| 



C = 



N K , 



£ trg ' A - C e ■ <J e 



(2) 



5. Efficiency 

A MC simulation is used for the determina- 
tion of the detection efficiency. In the e + e~ — > 
PPil) generator, corrections for initial state radi- 
ation [T3|, the Coulomb effect [T£][IB|, an d final 
state radiation |17| have been taken into account. 
The correction factors for these items, 1 + 8, F c , 
and Ff, respectively, are listed in Table □ 

For each energy point, 50,000 MC events are 
generated. MC events must satisfy the same se- 
lection criteria as used for the real data. The 
detection efficiencies, s, are given in Table 

6. Systematic errors 

Systematic errors come from uncertainties in 
the detection efficiency, trigger efficiency, lumi- 
nosity, and background contamination. The de- 
tection efficiency uncertainty includes the MC 
statistical error and the differences in the particle 
identification (PID) and tracking efficiencies for 
the MC sample and the real data. 

There are few e + e~ — > pp events, so p and p 
samples from J/tp — > ir + Tr~pp are used for the 
PID efficiency study. The momenta of the can- 
didate p and p tracks are required to be within 
30 MeV/c of the expected values for each en- 
ergy point. The samples arc obtained without 
using PID. The PID efficiency for each energy 
point is then detemined by the fraction of p and p 
tracks that pass PID selection criteria. The same 
method is used for MC J/tp — > Tr + ir~pp events 
to determine the PID efficiency for the MC data. 
At 2.0 GeV, there is a large difference in PID ef- 
ficiency between the MC sample and the data, so 
for all energy points detection efficiencies (F e ) are 
corrected for this difference, and the errors in the 
PID efficiency difference are taken as a source of 
systematic error. 

For e + e~ — > pp, possible backgrounds are from 
e + e~ — > e + e~(~f), /i + /i~(7), 7T + 7r~, K + K~ , and 
ppTT . MC events are generated for these five de- 
cays at s/s = 2.2, 2.6 and 3.0 GeV to estimate 
the amount of background contamination, which 
is included as a systematic error: 1.5% for 2.0 
and 2.2 GeV, 4.4% for 2.4, 2.5, and 2.6 GeV, and 
7.8% for other energy points. 

Systematic errors are listed in Tabled An un- 
certainty of 1.0% is taken for the initial state ra- 
diation correction. 
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Table 1 

Summary of results at center-of-mass energies from 2.0 to 3.07 GeV. N is the number of e+e~ — > pp 
events, C is the integrated luminosity, e is the detection efficiency, F e is the detection efficiency correction 
factor from the PID efficiency difference between the MC sample and the real data, 1 + 6 is the initial 
state radiation correction factor, F c is the correction factor for the Coulomb effect, Ff is the final state 
radiation correction factor, <7o is the measured lowest order cross section, and \G\ is the form factor. In 
the last two columns, the first error is statistical and the second is systematic. 
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45.8 ± 1.4 
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.02 


0.98 


64tg ± 4 


72±H ± 2 


2.5 


c+2.8 
— 2 2 
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.02 


0.98 


201±gi 3 ± 13 
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2.6 


24 ±4.9 
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.51 ±0. 
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54 ± 6 ± 2 
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0.96 


1 


.27 
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.01 
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27l^ ± 1 



Table 2 

The relative systematic error (%): Ae/e is the contribution from detection efficiency, including the MC 
statistical error (errAr), the PID difference (err p id) and tracking efficiency difference {err track) between 
MC and real data; Aetrig/strig is the contribution from the trigger efficiency; AC/C is the contribution 
from luminosity; and BG is the contribution from the background contamination. 



y/s Ae/e Aetrig/strig AC/C BG total 



(GeV) 
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errp/D 
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efficiency, and e tr i g is the trigger efficiency. The 
iowest order cross section is determined from 



IGI 



CTo 



(l + 6)-F c -F f ' 



(4) 



where 1 + 8, F c , and Ff are correction factors for 
initial state radiation, the Coulomb effect, and 
final state radiation, respectively. 

The form factor can be calculated from the the- 
oretical lowest order cross section jS] 



ff o = — (1 

3s 



2m 2 

^)|G( S )| 2 , 



(5) 



in which a is the fine structure constant, v is the 
proton velocity, m p is the proton mass, and \G\ 
is the form factor assuming \Ge\ — \Gm\- 

The trigger efficiency e tr i g for the hadron is 
about 100%, and the error is estimated to be 

0.5% hei. 

The cross section of e + e~ — > pp and proton 
form factor have been measured for 10 center- 
of-mass energies between 2.0 and 3.07 GeV. The 
measured values are listed in Table Q There is 
no signal found at ^Js — 2.9 GeV, but upper lim- 
its on the cross section and the form factor at 
the 90% C.L. are given, using the method from 
Ref. Ql|. 

For large momentum transfers pQCD pre- 
dicts [20] that q 4 G should be nearly proportional 
to the square of the running coupling constant for 
strong interactions, a 2 (q 2 ), yielding the relation 



\G\ = 



C 



s 2 ln 2 (s/K 2 



(0) 



where A = 0.3 GeV is the QCD scale parameter 
and C is a free parameter. In Fig. [3] BES re- 
sults are compared with other experimental pro- 
ton form factor results. The line is the |G(s)| 
energy dependence obtained by fitting all mea- 
surements with Eq. [B] and the result is consistent 
with the pQCD prediction. 
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Figure 3. Form factors measured by BES and 
other experiments. The line shows the energy de- 
pendence of | G(s) | by fitting all measurements. 
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